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Introduction 

The equivalent current concept has been applied implicitly in GTD 
solutions by various authors with Ryan and Rudduck^.g being among the 
earliest. It was formalized by Ryan and Peters [3,4] and applied to 
compute the radiation and scattered’fields from a number of radiating 
systems. It has since been applied to obtain the complete radiation 
pattern from reflector antennas[5,6]. 

Senior and Uslengi criticized the equivalent current concept in a 
paper discussing the scattering of a cone[7]. Burnside and PetersfB] 
corrected the misleading impressions generated by Senior et.al. and 
included the equivalent current in treating the second order diffrac- 
tion to improve the analysis of axial backscatter for a finite cone. 

The basis for the equivalent current concept is described compactly 
in our response to Senior's comments and it is repeated here. Fol- 
lowing that discussion some new previously unpublished results obtained 
using the equivalent current technique are presented that in essence 
fixes the minimum size conducting body such as a cone that can be 
treated using the ray optical format. 


Theoretical Background 

Consider a two-dimensional wedge whose edge is y directed and lies 
at the coordinate origin (r = 0). The diffracted fields for a plane 
wave at normal incidence is given by 
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where the parameters are illustrated in Fig. 1 with 

(J)± = ^ + ,1,0 . 


The incident and diffracted fields are given by 
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respectively. For the observation point removed from the shadow bound- 
aries and s sufficiently large 
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This form is used for all computations in this paper* Greater accuracy 
can be obtained for small values of s and for observation points much 
nearer the shadow boundary using a form given by Hutchins and 
l<buyoumjian[9] . Equation (2) can also be written in the form 
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It is seen that Equation (3) has the same range dependence as a line 
source. 


by 


The two dimensional fields of a y directed line current are given 
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for an electric line source with current , and 
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for a magnetic line source with current ^ as given by Harrington[10]. 
Equations (1) through (5) can be solved to give the equivalent currents |3]> 


(6) ^ 4*^ [G(n,<|)-) - G{n,ci>'-)] Ej 

and 
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Fig. 1 --Wedge Geometry 


If the shape of the edge is distorted, then the diffracted fields 
are predicted from these currents, positioned in space to conform to 
the shape of the edge, in that diffraction is a localized phenomena. 

The present study is restricted to the fields on the back surface 
or base of a cone (or disk) vjhen a plane wave is axially incident on 
the cone as shown in Fig. 2. The Geometrical Theory of Diffraction 
(GTD) (or ray optical analysis) would predict that these rays are dif- 
fracted across the back through center of the base to the opposite side 
of disk. A 90® phase jump is introduced after the ray passes through 
the center. Since all rays converge at the center the ray optical 
analysis would predict that the fields at this point become infinite, 
and this represents the caustic region. These rays are then diffracted 
at the opposite edge of the cone base and are designated as second 
order diffracted rays. The ray that is launched in the direction of 
source contributes to the back scattered field[ 8 ]. It is the fields 
that exist on the base of the cone that are the subject of the present 
study. The GTD analysis yields 

-E® = sin (J. ]jj- 2 V 0 ^n, a-r, 




-Cone Geometry 
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and 


for 0 < 4> < TT and 0 4*' < it 

or r < (j) < 2ir and ir < 4>' < 2tt 


-E® = -El 



a+r, 




for IT < 4>' < 2tt and 0 < <|) < tt 
or 0 <(])'< tt and ir < 4> < 2tt 


It is the factor that causes the GTD result to approach infinity as 

r 0. The sin 4> factor is introduced by the polarization of the 
incident wave. It would not be present in the acoustical case for the 
hard boundary condition. 

The equivalent currents have been introduced for the express pur- 
pose of obtaining the correct value of the fields in the vicinity of 
the caustic. 

The value of I®, given by Equation (6), required to compute fields 
on the surface of the conducting wedge is identically zero. Thus only 
is used in this paper since our interest is focussed primarily on 
the fields at the surface. Similarly for the acoustic case, only the 
hard boundary condition is considered. 

The field configuration 'of interest is that of the normal electric 
field in the vicinity of the caustic on the surface of a perfectly con- 
ducting body with a wedge type of discontinuity. A circular geometry 
is considered since it can be done analytically and the essential 
features of the caustic region can be established. This circular 
geometry could be used directly for studies of both cones and disks. 

For plane wave incidence along the axis of symmetry, the caustic v/ould 
exist in the vicinity of r = 0. Our attention is further restricted 
to the surface at the base of the cone corresponding to z = 0. The 
incident field is specified by 


P = 


y e 
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Analysis 

The equivalent current at the rim of a cone or disk (at r = a) is 
of interest and is given by 

— [G(n,(^ ) + G(n,(f) )J (-sin 
k 

j4 / 3ir \ 

(8) " iT 1”) 

for the geometry shown in Fig. 2. The z component of the electric field 
at the surface is given by 


(9a) 


jka 

4tt 
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It is most important to note that Equation (9) is valid only in the 
vicinity of the caustic. The current I'a is not an actual current but 
is only an equivalent current and is a valid representation only in the 
direction of the diffracted rays. Equation (9) will yield accurate 
results in the general vicinity of the caustic or for r small for our 
geometry. The equivalent currents are based on the concepts of edge 
diffraction and it represents a valid approach only in the vicinity of 
a converging ray bundle. It would yield approximate results for other 
values of r as long as a diffracted ray exists at that point, since this 
would represent a dominant term which would be obtained from Equation (9) 
via a stationary phase process. Balanis[13]has used the stationary 
phase concept to properly evaluate the fields caused by an edge diffrac- 
tion process. However, Equation (9a) cannot be used in general to 
evaluate a low level field at a point far removed from the caustic where 
there is no propagating ray. To the extent that Equation (9a) yields a 
result that deviates from the ray optical solution for an observation 
point far removed from the caustic, it is in error. The limits of inte- 
gration of Equation (9) give the total electric field on the base of the 
cone. This will set up a standing wave pattern since the currents for 
0< ())' < Ti will launch a set of diffracted fields as will the currents 
for ir< ((»' < 2¥. These two sets of diffracted fields will interact and 
create a standing wave pattern which will obscure the caustic field 
structure. If the integration limits are 0< i|i' < ir, only a single 
set of caustic fields are computed and no standing wave pattern appears. 
Most of the results presented in this paper are based on this range of 
integration so that the interference pattern is not included or 
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Approximate Form of Caustic Fields 

It is instructive to examine the characteristics of the caustic on 
a functional basis to obtain a further understanding of its form and 
extent and also to extract the phase jump introduced as a ray propagates 
through the caustic region. 

For observation points on the surface, the approximations 


(10) 

and 

( 11 ) 

( 12 ) 
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are made for r « a. Thus 
(13) 
where 
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By use of the relation 


gjkr cos((l>-4)') , ^ jn e3^<f>-V) ^ 

n=-oo 

-jk(a+|^') / \ 

(15) -Ee = Ez ^ 2j sin <}> J] (kr) e '' ®''G(n, 

The fields removed from the surface would be obtained using essentially 
the same procedure but v/ould require a contribution of the equivalent 
electric current of Equation (6). The approximation of Equation (12) 
is accurate within 10% provided r/a < 0.2 as is easily shown by means 
of the binomial expansion. Hov/ever, this may not prove adequate in the 
phase term and in fact for a/X = 10, r must be less than X in order 
that the phase error introduced by this approximation be less than 20°. 
On the other hand, for a/X =1.0 and r = 0.2 the maximum phase error is 
less than 10°, 
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Another basic limitation exists in that the equivalent current 
representation which is based on the asymptotic form of the diffraction 
coefficient requires ks (1 + cos (()-) » 1, This condition is well sat- 
isfied provided $ > X/2 except near the shadow boundaries = rr). 

This possibility does not occur here. 

The factor Oi(kr) has a maximum at kr = 1.8 and a minimum at 
kr = 3.8 corresponding to r/X = .27 and O.Gl, respectively. Later it 
is demonstrated that the two point diffraction mechanism is all that is 
required to properly evaluate the fields for r/X > 1/2 v/hich corresponds 
closely to the first minimum of 0](kr). For r/X > 1/2, the observation 
point then is outside the caustic region. 

Equation (15) could be obtained from the GTD analysis provided one 
uses the asymptotic form of 



The GTD solution has the form 
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which can now be related to Equation (15). 

Charge and Current Densities 

The magnetic field intensity at the suface of the conductor, is 
given by 


H = 
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jwy 


V X E 


ilormal 

Thus 


H is identically zero over the conducting surface, so Hot = 0. 

« 1 0=90° 
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It v;ould appear that we must rtow evaluate the additional components 
of the electric field intensity. However, these can be approximated 
rather accurately using the terms that have already been obtained. 

First, note that the equivalent electric currents need not be intro- 
duced since they would produce only a magnetic field perpendicular to 
the conducting surface at 0 = 90“. Second, and H<j, are related for 
a single ray by the intrinsic impedance of free space in the region 
where the ray optical computations are valid such that = + Zq H<J). 

This includes only a single ray propagating in one direction only and 
the sign is fixed by the direction of propagation. If there are two 
rays, is computed independently for each case and the total magnetic 
field intensity is, then, the vector sum of two terms. Finally in the 
vicinity of the caustic, under the condition that r « a, 

^ ^ 0 

Ep = E^ cos 8 

Eq ^ - E 2 sin 0 

where is given by Equation (13). With these approximations, the 
magnetic field intensity at the surface is given by 





Equations (15) and (16) may be used to evaluate the surface charge 
and current densities, respectively. The surface charge density is 
given by 


P = 


Co n 


• E 
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or 

y .2 

-jk(a + hj) . 3ttv 
p = 2e^j Oj(kr) sin (|>e G(^n, ^ J 

ThG surfacs currsnt dGnsity (0 n x H)is also rsadily obtained using 
Equation (16). 

Caustic Fields for Half Rim Excitation 

To obtain the field structure when only the rim for 0< < tt is 

excited. Equation (9b) was formulated. Now Equation (13) becomes 


E -E =Af s1n d*' 
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After some manipulation, this reduces to 
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The asymptotic form of the Bessel functions are 
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and recognizing the series[ll]as given by 
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plus the fact that 

(21 ) S((f)) = S(tt + (j)) 

for 0 < ()) < IT. These results are used to cast Equation (18) into the 
form 

(22) -Eg = A sin (j) e^('"""^) 

for 0 £ <fi £ TT 
and 

rr~ -irkr-^^') 

(23) -Ee = A \|:^ TT sin 4> e V 

for ir £ (}) £ 2 tt 

This last step v/ould apparently introduce an additional phase shift of 
-Tr/2 radians as the wave propagates through the caustic. However, the 
sin (() factor reverses sign and makes this +tt/2 radian phase shift which 
is in agreement with the previously accepted result. Note that Equa- 
tions (22) and (23) contain all of the restrictions implicit in Equation 


Straightforward application of ray optical techniques would give for 
distances removed from the caustic 


-Efi = G 


V 2 ; ✓ ^^kca-r) ' 


for 0 £ (j) £ TT 


and 
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for TT < 4> < 2fr . 

Comparing these equations with Equations (22) and (23) shows that the 
assumed phase jump has indeed been introduced after the ray passes 
through the caustic region. 
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Exp ressions for Fields at tlie Rim 

To date our attention has been focussed on the fields in the vicin- 
ity of the caustic. Let us now consider the fields at the rim of the 
disk. The excitation is again restricted to 0 < and the obser- 

vation point is now placed at r = a and tr < (}) < 2 tt. The range depen- 
dence is given by 


|r - r' 1 = va^ + r^ - 2ar cos ((})-<{)' ) 


It is again noted that Equation (9) is a valid representation only 
in the vicinity of the caustic. Any expression for the fields at the 
rim of the cone base can be considered to be valid only when the fields 
are large. Such an expression is formulated and evaluated at one point 
on the rim. It is used to indicate of the extent of the caustic fields. 
The integral has the form 


(24) 


-Efl = E, 


G^n » 2 


\ 

) / sin (|)' e 
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j2ka sm 


27T 
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sin 
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The integral remains finite at the edge because of the restrictions 
placed on <))’. Equation (9b) has b^n evaluated numerically for 
r = a, corresponding to observation points on the rim of the cone 
base. The results are to be discussed later. 

Equation (24) can be evaluated analytically for the specific value 
of (j> = TT. This point occurs at the position where ray optics would 
predict a field strength of 0. The solution of Equation (24) at (J> = i: 
is 



Equation (25) has its first null at ka = tt corresponding to a/X = 1/2 

This is considered to be the extent of the caustic region, i.e. when 

the base radius a/X < 1/2, the caustic covers the complete base of the 

cone, at least in the 4> = tt plane. At this stage, the postulate is 

made that the ray optics formulation is valid for a > X/2, Mote that 
■ Equation (25) is not accurate for a/X > 1/2. 
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Resul ts 


The results obtained using Equation (9a) and (9b) with a/X = 10 for 
X = 2 = 0 are compared with those obtained using conventional ray 
optical or GTD techniques in Fig. 3 and 4. Use of these 
equations eliminate the approximations made in the integration but those 
inherent in the equivalent current concept remain. The amplitude re- 
sults show reasonable agreement v/hen y is in excess of 0.2X for both 
the case where the entire rim of the base of the cone is illuminated 
(Fig, 3a) and also where only half the rim is illuminated (Fig. 3b). 

In these curves, the ray has passed through the caustic for negative 
values of y. Similar conclusions can be drawn for the phase of the 
fields in the vicinity of the caustic from Fig. 4. The effect of the 
propagation path has been removed in the phase plot of Fig. 5 to more 
dearly illustrate the caustic phase shift. This curve is offset by 
-45° before the caustic because of the e"J '"■/4 factor in the diffrac- 
tion coefficient of Equation (2). It ivould appear that the ray optical 
assumption of a 90° phase shift is essentially correct for y/X > 0.-5. 
This compares favorably to the value of a/X > 1/2 obtained in the 
preceding section pertaining to the base radius at which the caustic 
region covered the cone base. 

The conclusions thus have been made from results obtained using a 
cone with a 10X base radius. In the following paragraphs, the effects 
of the base radius are to be evaluated. 

Equation (18) indicates that the amplitude of the caustic is inde- 
pendent of the base radius provided the assumptions made in its deriva- 
tion are valid. Fig. 6 shows the amplitude curves of the normal E 
field for different cone base radii as obtained from Equation (9b). 

There is some dependence on the base radii but within the values of 
r,- a and X for which Equation (18) is valid this is a rather small 
variation, and indeed at r = 0, the magnitude (Ez| = .45 for all 
values of a. The black dots indicate the GTD solution. Note that 
the GTD solution is valid outside the regions illustrated in Fig. 6. 

The phase variations through the caustic region for different base 
radii are shown in Fig. 7, The 90° phase jump of the ray optical solu- 
tion is not achieved for smaller base radii. For a/X > 0.5, <j) = 180, 
a phase error less than 25° is made when the 90° phase jump is used. 
These results would be useful in the evaluation of charge density on 
the base of the cone or for computing the voltage at the terminals of 
an antenna element placed on the cone. 

In order to compute the back scattered fields of the cone, one must 
know the fields at the opposite rim of the cone. These have been com- 
puted using Equation (9b). As we have ascertained, the phase is the 
most significant parameter and it is shown in Fig. 8. The phase jump 
caused by passing through the caustic region is significantly different 
than the 90° phase jump of GTD for smaller radii bases. For a/X =0.5 
a phase error of 20° is introduced vdien the 90° phase jump is introduced 
according to the ray optical analysis. Fig. 6 would indicate that the 
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correct magnitude is predicted by GTD for r/A as small as 0.2. However, 
this indication seems to be valid only for observation points near 
({) = 270“. 

It has been noted earlier that the caustic introduces a significant 
perturbation fvom the fields computed by GTD at the rim when a/X < 1 / 2 . 
This has been verified again using Equation (9b) for a/X = 1 / 2 , 1, 2.5 
and 5.0 as a function of <|) when tt < i}) < 2 tt. As predicted previously 
the fields at cj) = tt are negligibly small. The GTD analysis would pre- 
dict the magnitude at the rim to be A( 4 >) = A(t}) = 270) sin <(). This 
assumption has been observed to be valid for 

no (|) (a/X = 0.25) 

210° < < 330° (a/X = 1/2) 

204° < (}) < 336° (a/X = 1.) 

200° < < 340° (a/X = 2.5) 

180° < <}) < 352° (a/X = 5.0) 

The phase term is usually the term in v^hich the error is first observed 
and the expected phase shift (that of Fig. 8 ) is obtained for 

252° < (fi < 288° (a/X = 1 / 2 ) 

210° < (J) < 330° a/X = 1 

206° < 4 . < 334° a/X = 2.5 

200° < <t> < 340° a/X =5.0 

It has been previously postulated that the extent of the caustic 
region is essentially r <. 0.5X. Yet it has been shown that the conven- 
tional ray optical solution produces a reasonable amplitude along the 
X = 0 axis for r >:0.2X. In order to resolve this question, the fields 
were computed atr=a, tr<(fi<2iT (on the rim) for a = X/4. It was 
not too surprising to find that the amplitude did not follow the 
A{(j)) = l\{4> =270°) sin ()> formula but were nearly constant ranging from 
0.235 at 4 ) =270° to 0.287 at 4) = 180 and 360°. This merely indicates 
that the observation point is nov; in the caustic region over the entire 
base as has been suggested previously even though the GTD analysis pre- 
dicts the correct amplitude at x = 0, y = a. Thus the fields are no 
longer local in nature and GTD does not give precise results for this 
small a cone base radius. 

Finally, contour plots of the magnitude of normal E at the base of 
the cone are shov/n in Figs, 9 and 10 as obtained from Equation (9a) and 
(9b), respectively. The standing wave pattern shown in Fig. 9 is caused 
by waves being launched at opposite sides of the base of the cone masks 
the shape of the caustic field. The caustic field region is more 
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clearly seen in Fig. 10 where only the currents for y > 0 are included. 
The ripple observed in y = 0 plane is in this case a valid result. 

Here the source distribution is at a radius of lOX and we are examining 
the fields in the vicinity of the caustic. It is also seen that the 
null position in this plot occurs at x = X/2 again showing agreement 
with the previous statement that the caustic occupies the region 
r/X < 0.5. 


Conclusions 


The caustic region on the back of a conducting cone for an axially 
incident electromagnetic wave is shown to extend to a distance on the 
order of X/2 from the center of the caustic. This fixes the applica- 
tion of GTD analyses for cone back scatter to ka > 3.0. Such compu- 
tations[8] have been made for cones as small as ka = 2 and reasonable 
results have been obtained. It may be possible to improve these results 
by replacing the 90° phase shift commonly caused by the caustic by the 
value given in Fig. 8 plus 45°. However, smaller cones could not be 
treated since the caustic v/ould cover the entire base and the fields 
no longer have the local nature required by the ray optical or GTD 
analyses. 

As has been noted earlier the results of this paper are applicable 
to finding the response of elemental antennas placed on the cone base, 
even for antennas mounted at the center of the caustic region. It is 
possible to consider smaller cones than the back scatter analysis, 
probably for base radii as small as X/4, as fixed by the restriction on 
this distance of the observation point from the edge of the cone rim. 
Second order diffractions would again be in error but the fields from 
these rays would be small compared to the first order diffracted fields. 

Similar comments may be made concerning the surface charge density 
induced on the cone base by an incident wave. 


Appendix. Equations for 'Acoustic Caustic" 
using Hard Boundary Condition 

The pertinent caustic equations are presented here for the hard 
acoustic boundary condition. 

Equation (8) becomes 
(A-,) ih.aG(n,^) 

Equation (9a) becomes 
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Equation (13) becomes 

u = A f d(j)' 

0 

Equation (15) becomes 



Equation (A-4) differs significantly from Equation (15). The polar- 
ization factor (sin <!)') in Equation (13) caused the coefficient of the 
jQ(kr) term to be zero for the electromagnetic case. If we retain the 
concept that the caustic region extends to the first null, then the 
"acoustic caustic" for the hard boundary condition extends only to 
r < 0.4X as compared to r < 0.61 for the electromagnetic caustic. 

The form obtained for the restricted ranges of integration equivalent to 
Eq. (9b) is of limited merit in this case because of the severe discon- 
tinuity in the current density at c|)' = 0, tt. This difficulty is averted 
in the electromagnetic case because of the presence of the sin 4>' factor 
in the integrand of Equation (9a). 

It is worth noting that by use of the asymptotic form of On(kr), 
both Equation (15) and Equation (A-4) can be obtained from the GTD 
analysis. To obtain Equation (A-4) for example, the expression for the 
GTD analysis is 


u = u’ ^ ' 2v0(n, a-r, 2vB(n, a+r, e"^ 2 

and using the far field form of the diffraction coefficient yields 

/ 37r\ r . . it . jr 

^ /2irka 


i -okaf 

= u* e N 


/2irkr ^ 



= U’ 2G(n, f ) g cos(kr - J) 


where 


l/llr ■ f) 
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